home *** CD-ROM | disk | FTP | other *** search
/ Pascal Super Library / Pascal Super Library (CW International)(1997).bin / MATH / NRPAS13 / SIMP1.PAS < prev    next >
Pascal/Delphi Source File  |  1991-04-29  |  740b  |  29 lines
  1. PROCEDURE simp1(a: glmpbynp; mp,np,mm: integer;
  2.       ll: glnparray; nll,iabf: integer;
  3.       VAR kp: integer; VAR bmax: real);
  4. (* Programs using routine SIMP1 must define the types
  5. TYPE
  6.    glmpbynp = ARRAY [1..mp,1..np] OF real;
  7.    glnparray = ARRAY [1..np] OF integer;
  8. in the main routine. *)
  9. LABEL 99;
  10. VAR
  11.    k: integer;
  12.    test: real;
  13. BEGIN
  14.    kp := ll[1];
  15.    bmax := a[mm+1,kp+1];
  16.    IF (nll < 2) THEN GOTO 99;
  17.    FOR k := 2 TO nll DO BEGIN
  18.       IF (iabf = 0) THEN BEGIN
  19.          test := a[mm+1,ll[k]+1]-bmax
  20.       END ELSE BEGIN
  21.          test := abs(a[mm+1,ll[k]+1])-abs(bmax)
  22.       END;
  23.       IF (test > 0.0) THEN BEGIN
  24.          bmax := a[mm+1,ll[k]+1];
  25.          kp := ll[k]
  26.       END
  27.    END;
  28. 99:   END;
  29.